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ABSTRACT 


It is important for a student of Space Science to have 
the opportunity to thoroughly understand the principle of an 
artificial satellite orbit. This thesis consists of four 
computer graphics programs that will enable the student to 
see what an orbit is and how it works. The first program 
demonstrates the shape of an orbit in two dimensions 
resulting from initial altitude, speed, and flight path 
angle. The second program draws an orbit in three dimensions 
around a sphere based on the input of the classical orbital 
elements. The third program traces the ground track of a 
satellite over a map of the earth. And the fourth plots the 
ground track of a geosynchronous satellite over a map of the 
earth. The student can learn about orbits by entering the 


orbital elements and viewing the resultant orbit. 
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I. INTRODUCTION 


In Principia Mathematica Sir Isaac Newton theorized 
Pte by firing a Cannonball from the top of a high mountain 
at increasing speeds it would travel farther and farther 
Geer! at last it esnmand make a complete orbit around the 
earth [Ref. l:p. 66]. As with many great discoveries, the 
idea of orbital flight came from a simple idea. Three 
Tetieed years later mmonmsm0ctober 1957, the launch of the 
first artificial satellite around the earth demonstrated 
that Newton's theory was essentially correct. 

Spaceflight has always been tied closely to computers. 
Graphics is defined by the Oxford English Dictionary as "of 
Cadaver os = pDalntine, engraving, .etching,etc.; vividly 
descriptive, lifelike; of diagrams and symbolic curves." 
Traditionally computer graphics evolved as a means of 
interpreting and displaying the increasing volume of numbers 
which the modern computer can manipulate with tremendous 
rapidity and accuracy. In 1963 Ivan Sutherland developed the 
first interactive computer graphics system called the 
"Sketchpad". 

This advance in computer graphics enabled a user to interact 
directly with a computer graphics program in real time. 


Ret. 2p. vii | 


A. OBIECl int 
The laws of celestial mechanics) as tormulated@)y sl yee 
Brahe, Kepler, and Sir Isaac Newton came from simple ideas 
and meticulous observations of the heavens above. The 
intent of this thesis is to use the graphics capabilities of 
the microcomputer to help students of space science 
understand the simple laws of orbital mechanics that pertain 
to artificial satellites. It is sometimes easier to 
understand a scientific principle by first observing a 
demonstration of what it does and then studying all the 
theories, equations, formulas, and laws that explain the 
principle than by trying to visualize what it does fron 
studying the equations alone. 
Four software programs are the result of this effort: 
1. INITIAL CONDITIONS demonstrates the shape of the orbit 
in two dimensions after the user enters the initial 
altitude, the initial flight path angle, and either 
the initial speed or the desired apogee height of the 
satellite. 
2. ORBIT plots an orbit around a sphere representing the 
earth after the user inputs the classical orbital 


elements. 


3. GROUND TRACK traces the ground track of an orbit ona 
modified cartesian projection map of the earth. 


4. GEOSYNC displays the ground track of a geosynchronous 
orbit of any inclination or longitude of the equator 
crossing on a modified cartesian projection map of the 
earth. 


Be. HARDWARE 

ire Weeislon FOmiscuagmticrocomputer for this project 
was based on the idea that microcomputers are more 
accessible and easier for the average student to use than 
ane either neta yonivprinelss or mainframes. Ung i addres) y 
recently the practice of computer graphics was restricted to 
special research laboratories and large scientific computing 
facilities. This was because the equipment necessary for 
computer graphics was expensive and bulky in size. In the 
mid-1970's came the flashing lights of the video arcades and 
inexpensive home video games. Several years later computer 
graphics became more available as the microcomputer 
revolution erupted. Apple computer offered the first 
microcomputer that could be operated with a minimum of fuss. 
And it had the added advantage of several types of color 
eraphics. TRS-80, PET, Northstar, Heathkit, and in the late 
1970's, Atari and Texas Instrument improved on what the 
Apple could do with graphics. In the summer of 1981 IBM 
entered the competition with the IBM PC. With Microsoft 
supplying the software IBM offered a microcomputer that 
quickly carved out a large share of the market. [Ref. 3: 
pp. 6-8] 

In order for a computer to be useful as a learning 
tool, it must be readily available and easy to operate. The 
IBM Personal Computer (PC) was selected for this project for 
those reasons. It is one of the more popular PC's and is 


widely emulated by other computer manufacturers. A major 


portion of the programming for this project was done ona 
Compaq computer, a popular IBM compatible. 

A typical microcomputer OT PC IsM@ontig tired sae 
central processing unit (CPU) which is the "brain" of the 
computer, a keyboard tee is the means to input user 
intructions to the CPU, and a screen (CRT) which provides 
the output to the user. In order to support computer 
graphics the PC must have a color/graphics adapter that 
produces both the color and the graphics on the screen Jae 
store the programs and the software used to develop the 
programs a disk drive or other memory storage device is 
necessary. It is preferable to have two disk drives to cut 
down on the number of times a floppy disk has to be inserted 
and removed. The last recommended hardware component is a 
printer, It is very helpful to print out a program in order 
to debug it. And with the proper graphics software and 
printer it is possible to print the graphics output from the 


screen to the printer. 


Cy SOFTWARE 

The software that controls the operation of a computer 
is called the operating system. It directs such operations 
as loading software to the internal memory and saving 
programs to the floppy disks. PC-DOS is one of the 
operating systems available for the IBM PC. It was chosen 
for this project because is offers the computer language 


BASICA. IBM BASICA was developed by Microsoft, Inc. and was 


derived from BASIC which was developed at Dartmouth College 
to be a computer language that is easy to learn and easy to 
work with. It is normally an interpreted language and that 
means that it is edited and programs are run in the same 
working environment. This makes BASIC a very user-friendly 
computer language. The disadvantage is that it is 
relatively slow and graphics can take a lot of computational 
fame betore being@displayed. Compilers are available to 
speed things up by taking the source code and Priniacel a best 1 
a form the computer can handle and display much faster. 
IBM's BASICA offers some of the best graphics commands 
available on a microcomputer. The Graphics Definition 
momeuace (GDL) alone with the LINE, CIRCLE, PSET, and PRESET 


commands provide powerful tools for graphics programming. 


is PROGRAM DESIGN 

The general form of any software program is to accept 
user input from the keyboard, make the necessary 
calculations, and then output the results. In designing 
these programs an attempt was made to make them as user- 
friendly as possible. To start the programs, the user 
inserts the disk in the "A" or default drive and turns on 
the computer. The INTRO program is automatically loaded and 
run by a program called AUTOEXEC.BAT. Lema aLrrerent 
version of DOS is used, it is necessary to make sure that 


the BASICA is version 2.1 or higher. It might then be 


necessary to type BASICA at the system prompt and then 
‘load"INTRO",r’ to start™@the progransc. 

BASICA has several features to help set up a graphics 
screen. The medium resolution screen is 320 by 200 picture 
elements (pixels). In additions fo va aera color there 
are three colors that can be used for drawing on the medium 
resolution screen. A higher resolution screen is available 
but it will not support color graphics. A partilongo tence 
screen can be used as a “window” by using the WINDOW and 
VIEWPORT commands. This allows many graphics functions such 
as scaling the size of objects, clipping objects Ehatuare 
too big for the window, and clearing only the screen inside 
the window. 

Once the screen is set up, the user input is entered 
from the keyboard in response to prompts from the computer 
program. The entry is usually a number or a one or two 
letter response. The number can be entered in either real, 
integer, or floating point format. A yes or no response can 
be entered by typing a "y" or a "n", To select a program 
from the menu, type the first two letters in the name. En 
most cases the program will doalimited check to see if the 
input is appropriate. If not, a beep will sound and the 
question is again asked. 

After the initial input is accepted, the program will 
do the calculations for the orbital data. The program then 
enters a loop that plots the orbit from the perigee point 


through one complete orbit around the earth and back to the 


12 


famicee pointe. lhe Orbitwassplotted in steps of 2 time pi 
divided by an integer, i.e., 2 times pi divided by 50. 
After the loop is completed, the program will offer the user 
a choice of drawing another orbit or returning to the menu. 
Depending on the See iin the orbital data is either given 
as the orbit is being plotted or as a viewing option after 
the orbit is drawn. 

The physics formulas and equations in the programs came 
from a variety of sources. The BASICA language allows the 
programmer to enter an equation in the program in pretty 
much the same form as it is written in a textbook. Several 
exceptions to standard notation include using "*" for the 
multiplication operation, sqr( ) for the square root, and 
expressing exponentials as "X*2" for X to the second power. 
The equations in this thesis are written in the BASIC 
algorithmic language. The reason for doing so is to make 
the program listings as easy as possible to understand. The 
equations in the descriptions of the programs might not 
exactly match the equations in the listings because of 
special programming techniques which must be included in 
most computer programs to handle such problems as "division 
by zero". Another reason is that only a few trigonometric 
functions are directly computed in BASIC and the rest must 


be derived from those. 


There are several common assumptions and constants used 
throughout the programs: 

1. Both the earth and the satellite are considered to be 
spherically Symmetric. This allows these bodies to be 
treated as though their masses are concentrated at 
their centers (point masses). 

2. There are no other forces acting on the bodies other 
than the gravitational forces which act along a line 
joining the centers of the earth and the satellite. 

3. Equatorial radius of the earth = 6378.165 kilometers 

4, Gravitational constant of the earth = k = GMe = 
398603.2 km?/sec* [Ref. 4:pp. 11-12, Ref. 5:p. 113, 


Ret..,67 0. 374.) 


II. PROGRAM INTRO 


This program displays the title page, introduction, 
and lists the menu of programs. It also ties all the 
programs together so that they can be called and run from 
the menu. Program INTRO is automatically loaded and run 
when the disk is inserted in the "A" or default drive and 
the computer is turned on. I[f the computer is already on, 
the system can be rebooted by pressing the Cntr, the Alt, 
and the Del keys at the same time, and releasing them 
together. 

After the title page and a brief explanation of the 
programs, the user can run any of the programs by typing in 
the first two letters of the program name. The INTRO 
program will then load and run the chosen progran. After 
the program is completed, the user has the option of running 
it again or returning to the menu for another choice of 
programs. 

When the user wants to terminate using the programs, 
the only precaution is to make sure that the red light on 
bhe disk aqrive is out, indicating that the drive is not in 
use. The disk can then be taken out of the disk drive or 


the computer can be turned off. 


III. PROGRAM INITIAL CONDITIONS 


This program demonstrates the physics of the two-body 
problem. A similar program was written by Dr. John N. Dyer, 
Distinguished Professor eor ss woeec, Named Postgraduate 
School, for an early model Hewlett-Packard computer. The 
idea for this program and the majority of the equations come 
from Dr. Dyer's program; any equations used that were not 
ErLOm Sthivs) SO mec will be referenced. After an 
introduction, the user is prompted to enter the altitude 
above the surface of the earth of the satellite in 
kilometers. The altitude of the satellite is added to the 
radius of the earth to get the initial altitude of the 
satellite (Ro). The program then asks for the initial flight 
path angle in degrees (Bo), which is the angle between the 
velocity vector and the local horizontal plane. The final 
user input for the program is the selection of the speed in 
kilometers per second (Vo), or the desired apogee altitude 
above the surface of the earth in kilometers. If the user 
selects the apogee altitude, the program calculates the 


initial speed (Vo) as follows: 


Ree oy Rene 


XXX = ((2-X-X*2*cos(Bo)*2-sqr((2-X-X*2*cos(Bo)*2)*%2-4 
*(1-X)*(1-X*2*cos(Bo)*2)))/(1-X*2*cos(Bo)%*2) 


Vo = sqr((XXX*K)/Ro) 
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Figure 1. 


The program uses the initial conditions <oetimaeeme 


2 


specific angular momentum (SAM) in km“ per second 


SAM = Ro * Vo * cos(Bo) [Ref. 8:p. 74] 


The specific mechanical energy (SME) in km? per sec* is 


SME = (Vo*2/2)-(k/Ro) [Ref. 8:p. 74] 


The shape of the orbit is determined Sy the eceentricizs 


The equation for eccentricity is 


AAA = (Ro*Vo0“%2)/k 
ECC = sqr((AAA-1)*2*(cos(Bo))*2+(sin(Bo))*2) 
Tf the eccentricity = 0; the ordlt lswancenct ec. 


If the eccentricity is between O and 1, the orbit is an 
ellipse. 
If the eccentricity = 1, the orbit is a parabola. 
If the eccentricity is greater than 1, the orbit is a 
hyperbola. 
If the eccentricity is less than Q, the orbit is an ellipse 
with the apogee at the initial condition point. 

Two velocities of interest are calculated next. They 
are the circular orbit speed (VCIR), and the escape speed 


(VESC) in km/sec 


VCIR sqr(k/Ro) 


VESC 


sqr((2*k)/Ro) [Ref. 4:pp. 34-35] 


-— 4 


at « ea 
i a PS LP ed j % 
ky mee rs A v 
: f 4X, | | ‘, 
au j { 
‘ ie Ms TE ro ABE. PL a a fs ioe of * age ’ po ae 
x i 
i . : 
ee es 
- oh : = = 
44 on — 2 958 tet. 6 OF . “2% 
ai wees SFG cK PU Ct 
an cent ee : =: 
3: Binet ee Oe ees 
Pe nate AS n- cca ul te ee e 
ie Sere ees oa Pa 


Se i 8 came fk OM OO ee ee eB ee ee ee i 0 i ee _ ee. ee: ee ON ee ee et ee oe EE. eee __J 


Figure 2, Geometrical dimensions of an orbit 


If the orbit is an ellipse, the following orbital data 


is also calculated: 


The semi-major axis(A), the semi-minor axis(B), and the 
semi-latus rectum (PP) in km 


A = abs(Ro/(AAA=-2)) 
B = sqr(A°2*(1-ECC“%2)) 
PP= A*(1-ECC*2) [Ref. 7:pp. 97-98] 


The radius at apogee (RAP), and the radius at perigee 
(RPE) in ko 


RAP = A*¥(1+ECC) 


RPE = A*(1-ECC) 


The velocities in km/sec at apogee (VAP), and at 
perigee ( VPE) in km/sec 


VAP SAM/RAP 


ll 


HT 


VPE SAM/RPE [Ref. 4:p. 27] 
The period (1 ){ 820, minwees 


T = SOR((4* pi’ 2*A°3)/k) [Ret 3 4ep oo 
If the orbit is circlular, the folewine “data 


calculated: 


The eccentricity = QO. 
The semi-major axis is the initial radius (Ro = A = B = 
PEE 


The specific angular momentum (SAM), the specific 
mechanical energy (SME), and the period (T), using the 
same formulas as for an ellipse since the circle is 
special form of the ellipse. 


If the orbit. ais parabolic, themitcltovane. dataues 


calculated: 


Ihe eCcceneriecvteyn ECG nel. 

The specific mechanical energy (SME) = O. 

The specific angular momentum (SAM), SAM >D O. 
The altitude at senna (Ro-Re). 

The semi-latus rectum (PP). 


The escape speed (VESC). 


tite OlLbDL&e ts Myperbolic, the following is 


calculated: 


ime eccent paicastye( ECC). E@C.> 1. 

The altitude at perigee (Ro-Re). 

The specific mechanical energy (SME). 

The specific angular momentum (SAM). 

The escape speed (VESC). 

The hyperbolic excess speed (HES) in km/sec 


HES = sqr(Vo*2-VESC*2) [Ref. 4:p. 40] 


If the orbit intersects the earth's surface it is 
considered a ballistic trajectory and the following is 


calculated: 


The maximum altitude above the surface in km 
MAXALT = RAP-Re 

The true anomaly at launch in degrees 

THETAO = arccos((1/ECC)*(AAA*cos(Bo)*2-1)) 
The true anomaly at impact in degrees 


THETAI = arccos((1/ECC)*((AAA*Ro) /Re)*cos(Bo)*2-1) 


The range in km 
RANGE = Re * (2*pi-THETAO-THETAT ) 
' The eccentric anomaly at launch in degrees 
EO = arccos((1/ECC)*(1-(Ro/A)) 
The eccentric anomaly at impact in degrees 


EL = 2*pi-arccos((1/ECC)*(1-(Re/A)) 
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Ballistic trajectory 


Figure 3. 


Time of flight in minutes 

FLTTIME = sqr(A*3/k)*(EI-ECC*¥sin(EI)-EO+ECC*¥sin(EI)) 
Optimum launch angle in degrees for maximum range 

Bum 2/(1/(2*A-Re)+1/(2*A-Ro) ) 


BOPT = arccos(sqr(k*P1)/(Ro*Vo)) 


If the launch angle is 90 degrees, the trajectory is 
@omsidered a vertical launch and the following is 


calculated: 


The maximum altitude in km 
RMAX = (2¥*Ro)/(2-AAA) 


MAXALT = RMAX - Re 


The rise time in minutes 
Q = arcsin(sqr(Ro/RMAX)) 
TRISE = .5*sqr(RMAX*3/(2*k) )*( pi+sin(2*Q)-2¥*Q) 

If AAA is greater than or equal to 2, then the satellite 
escapes. 

After the orbital data is calculated, the program plots 
the orbit. A circle representing the earth is drawn using 
the CIRCLE command. The size of the circle is determined by 
Poe radwuseot the orb¥t (Ro). The Screen coordinates (X,Y) 
are determined from the spherical coordinates (R, THETA). 
A loop is used to plot the orbit 


FOR THETA = 0 to 2*pi STEP pi/30 


The radius (R) in spherical coordinates is determined 


Tf thevorbit 2s 
ealiete Sidlala © 

R = A 
elliptical: 


R = (A *® (1=ECC*2))/(1+4ECC*cos(iinis)) 
parabolic: 


R = (2*RPE)/(1l+cos(THETA) ) 

hyperbolic: 

R = (A *® (ECC*%2-1))/(1+ECC*cos( THETA) ) 
vertical launch: 


R = (2*Ro)/(2-AAA) where THETA = O. 


The screen coordinates, (X,Y), are calculated from the 


spherical coordinates 


X 


R * cos(THETA-THETAO) 


x 


R * sin(THETA-THETAO) 
These points are then plotted on the screen with the BASICA 
PSET command. 

At the end of the loop, the command NEXT THETA causes 
the program to return to the start of thesloop anti) jege 
orbit is completed. The user is offered the opportunity to 
plot another orbit over the first. It must be an orbit with 
the same initial altitude because the size of circle 
representing the earth is determined by the altitude of the 


orbit. At the end of every orbit the user can elect to view 


the data associated with that orbit. BEeeiiiVsmOoptLon 1s 
chosen, then another orbit cannot be plotted over the 


previous orbit. 


IV. PROGRAM ORBIT 


The program INITIAL CONDITIONS plots a two- 
dimensional orbit in which size and shape of the orbit are 
determined. Program ORBIT takes this one step further and 
shows the orientation of the orbit in three-dimensions. 
Five independent "orbital elements" are traditionally used 
to describe the size, shape and orientation of an orbit. A 
sixth element, time, is used to pinpoint the position of the 
satellite along the orbit at a particular time. The six 
classical orbital elements are: 


1. Semi-major axis (a) - a constant defining the sizeof 
the orbit. It is measured as half the distance from 
the perigee point to the apogee point. 


2. Eccentricity (e) - a constant defining the shape of 
the orbit. It varies from zero for a circle, between 
zero and one for an ellipse, one for a parabola, and 
greater than one for a hyperbola. 


3. Inclination (i) - the angle the orbital plane makes 
with the equatorial plane measured from the equator 
up in a counter-clockwise direction. 


4. Longitude of the ascending node (Q) - the point where 
the satellite crosses the equator in a northernly 
direction (ascending node) measured counter-clockwise 
from the prime meridian (0 degree point). 


5. Argument of perigee (w) - the angle, measured along 
the orbital plane in the direction of the satellite's 
motion, from the ascending node to the perigee point. 


6. Time of perigee passing (T) - the time the satellite 
is at the perigee point. The satellite position can 
be pinpointed by measuring the time elapsed from the 
time it passes the perigee point. [Ref. 4:p. 58] 
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Figure 4. 


The equations in this program came primarily Sirem 
Methods of Orbit Determination by P. R. Escobal [Ref. 6] 

Any exceptions will be referenced. 

ORBIT begins with a brief introduction followed by the 
user input of the first five classical orbital elements. Tf 
the eccentricity is zero, the program skips the input of the 
argument of perigee and longitude of the ascending node as 
they are not defined in the case of a circular orbit. In 
order to keep the program simple, it was decided to show one 
complete orbit from perigee point to perigee point. fhe 
sixth classical orbital element, time of perigee passage, 
would be measured from the perigee point. 

Once the semi-major axis is entered, the program sizes 
the sphere representing the earth and draws it in the 
viewport. 


The program converts the input from degrees to radians 


TORAD = .01745329 


X = X * TORAD 
distances are converted to canonical units 


Re 6378.165 km 


I 


ER A / Re 


i! 


ite Texte SECM an  Ehe program is to convert the 
coordinates of the satellite's position to the P-Q-R 


coordinate system 


W = argument of perigee in radians 
I = inclination in radians 
O = ascending node in radians 


oe Gos (W)*eCOS(O)-SINCW)*SINCO)*COSCE ) 
Bye GO5(W)*SIN(O)+SINCW)*COS(0)*COS(T) 


PZ = SIN(W)*SIN(I) 


moe —oINCW)*COS(0)-COS(W)*SINCO)*COS(T) 
QY = -SINCW)*SIN(O)+COS(W) *COS(0)*COS(T) 


mee= COS(W)*SINCL) 

the mean motion is calculated 
N= k / ER*1.5 

the period is 
Pee EWOPRL / N 


The main program loop divides the orbit into time steps 
of the period divided by 50. The loop goes from the perigee 


point for one complete orbit back to the perigee point 
HOR it =O 10 P STEP 50 


iiemrenastmcalemmuation within the loop is the solution to 


Kepler's equation. This is accomplished by a subroutine 


that uses another loop to find the eccentric anomaly (EN1). 
Lines 1620-1760 in the program listing contain the 
subroutine for the solution to Kepler’s equation. This 
method came from Dr. Fuhs' class notes eee 9]. 

The next step in ne loop is to convert the coordmnares 


to the XW-YW-RW coordinate system 


XW = ER*¥(cos(EN1)-ECC) 
YW = ER¥sqr(1-ECC*2)*sin(EN1) 
RW = ER*(1-ECC*cos(EN1)) 


the velocity is found 


VEL = sqr(k * ((2/(RW*¥Re))-(1/A))) [Ref. 9] 


the revolutions the satellite makes around the earth each 


day (N1) is 


Nl = 16.997 4°(1/ER) 125° [| Rete eee 


The program then prints the velocity in km/sec, the 
revolutions per day (N), the orbital period in minutes, and 
the radius of the orbit in km, on the top portion of Bene 
screen. 

If the orbit strikes the earthe(suborbital) eee 


program flashes a warning. 
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Example of Orbit Program 


Figure 5. 


The next step in the loop converts the coordinates to 


the X-Y-Z coordinate system 


X = XW*¥PX+YW*QX 

Y = XW*¥PY+YW*QY 

Z = XW*¥PZ+YW*QZ 
Before plotting the coordinates to the screen, a hidden 
line removal algorithm is used to determine those 


coordinates hidden behind the globe. They are not plotted 
and the loop continues, 

The CIRCLE command is used to plot the orbit around the 
globe. When the loop is finished the NEXT TTT command sends 
the program to the top of the loop until the orbit is 
complete. After the orbit is drawn, the option Of )draywasme 
another orbit over the first is given. The user can also 


elect to?’ rettrn to the menu. 


V. PROGRAM GROUND TRACK 


This program takes for user input the same set of 
classical orbital elements aS program ORBIT. It displays 
m@emeround trace of the satellite rer the surface of the 
earth on a modified cartesian projection map of the world. 
The map was drawn pixel-by-pixel and the entire screen saved 
with the BSAVE command. At the start of the program it is 
loaded to a non-viewed screen with the BLOAD command. The 
non-viewed screen is then paged to the viewing screen. 

This program is identical to ORBIT in converting the 
initial input to the X-Y-Z coordinate system. Instead of 
drawing a sphere, the program draws a map of the earth. The 
equations again are from Escobal [Ref. 6]. The exceptions 
are referenced. 


To get the longitude from the X-Y-Z coordinates 


Ol = ASCENDING NODE 
ROTEAR = ROTATION OF THE EARTH = .004369 rad/min 
ALPHA = arctan(Y/X) 


LONGITUDE = ALP -(CROTEAR*TTT)+01 


The most difficult programming problem of this project 
was to keep the longitude in the correct quadrant after 
taking the arctangent of Y divided by X. The solution to 
the problem is found in the subroutine in lines 1330-1410 in 


the program listing. 


The latitude is found from the X-Y-Z coordinate system 


R = sqr(X*2+Y°2+Z*2) 
LAlwomarecs ince waa) 

The program Changes tie mlaritude and longitude from 
radians to degrees and there are checks in the program te 
keep the longitude between O and 360 degrees and the 
latitude between +90 and -90 degrees. 

The latitude and longitude are then scaled to the 
screen and plotted from the perigee point using the CIRCLE 
command. GROUND TRACK has the option of drawing more than 
one orbit in order to see what the intertrace spacing looks 
like. The intertrace spacing is the distance on the earth 
between the first and subsequent ground tracks. 

After the ground track for a particular satellite is 
finished, the user can enter the orbital elements to plot 
the ground track of a different satellite or return to the 


menu. 
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VI. PROGRAM GEOSYNC 


A geosynchronous orbIlE is a Special §Cace Vo em 
circular orbit. The altitude must be-‘constant at about 
19,300 nm., 5.6 earth radii, or 35,750 km. The period is 23 
hours and 56 canteen The only variables are the 
inclination and the longitude of the equator crossing. If a 
geosynchronous satellite has an inclination of zero degrees, 
then it will appear to hover motionless over a point on the 
equator. “slioaws then called a geostationary satellite. The 
general pattern of the ground track of a geosynchronous 
satellite is a figure eight about the equator. 

Program GROUND TRACK can plot a geosynchronous 
orbit, however it is very slow in doing so. Dr. Allen E. 
Fuhs, Distinguished Professor of Aeronautics at the Naval 
Postgraduate School, presented an alternate set of equations 
to plot the ground track of a geosynchronous orbit [Ref. 9]. 
The computation time is much shorter, therefore the computer 
can plot it much quicker. 

GEOSYNC uses the same map as GROUND TRACK. After the 
user inputs the inclination and longitudinal crossing point, 


the program calculates the rotation rate of the earth 


We = (2¥*pi) / 24 radians per hour 
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the rotation rate of the satellite equals the rotation rate 


of the earth 


Wo = We rad/hr 


The time loop is in hours and goes from 0 to 18 with 


steps of .25 hours 

For T/ss0710 135i rae) 
The latitude is calculated 

LAT = ARCSINCSINCI)*SINCWo*T) ) 
and the longitude is 

LONG = ATNC(COS(L)*TANCWo*#T) )-We*T 


The latitude and longitude are then converted back to 
degrees, scaled to the screen, and plotted to the map of the 
world with the CIRCLE command. 

An additional factor that needed to be dealt with in 
this program was the adjustment necessary when the orbit was 
located in the 4th quadrant (X and Y are both negative). In 
this case another loop was necessary. The only difference 
between the main loop and the 4th quadrant loop was the 


addition of 2 times pi to the Yoneaitrtude. 


VII. RECOMMENDATIONS FOR FURTHER DEVELOPMENT 


These programs are meant to be used: by beginning space 
science students. They were kept at an-introductory level. 
As such, they should be considered as an initial attempt and 
improved upon. The next step from a simple Keplerian orbit 
would be to add the various perturbations that act ona real 
artificial satellite. These include the effect of other 
attracting bodies (the moon), atmospheric drag and lift, the 
oblateness of the earth, magnetic torque, at eerolian: 
radiation pressure [Ref. 4:p. 386] Perturbations must be 
considered to achieve reasonable accuracy as far as 
predicting satellite position over a long period of time, 

Another idea for further development would be to write 
programs for basic orbital maneuvers or rendezvous problems, 
muna , interplanetary, and missile trajectories could also 
be programmed. The equations for these programs are found 
in several texts on astrodynamics [Ref. 4, 6, 7, 8]. 

The Air Force Academy and The Naval Academy have 
orbital programs that run on mini-computers that support 
graphics. A good project would be to develop a similar 
capability at The Naval Postgraduate School. The graphics 
capabilities of the larger computers are much greater than 
those of most microcomputers. Two new microcomputers that 


have tremendous graphics capabilities are the Amiga and the 


Atari S$T-1040. Orbital eraphics programs written for thease 
two machines would be able to achieve better graphics than 
on most older microcomputers. 

A final recommendation is that the GROUND TRACK program 
could be modified to show ground coverage, number of 
satellites in a constellation, and other elements necessary 


for planning a real-world artificial satellite application. 


APPENDIX A 


INTRO PROGRAM LISTING 


may’ INTRO 

20 

30 'This is the introduction and menu program 
40 'that controls the other programs. 
or 

60 

70 'set up and draw the box around the screen 
so ' 

Pemeclo?SCREEN 1,0:KEY OFF 

100 COLOR 9 

ipo ’ 

120 'title page 

30 ' 

140 GOSUB 870 

iO ' 

160 ‘description of programs 

| 70 

180 GOSUB 610 

190 ' 

Z00 ' 

210 'draw the box around menu 
P20) LINEC(O,0)-(319,199),1,B 

230 LINE(10,10)-(309,189),1,B 
240 PAINT(5,5),1 

Zo 

260 'print the choices on the screen 
270 ' 

ASO LOCATE 5,17 

290 PRINT"PROGRAMS: 

BOO LOCATE 9,12 

310 PRINT"INITIAL CONDITIONS 

Peo eGUGATE 11,12 

Seo PRINT" ORBIT 

eeeo LOCATE 13,12 

350 PRINT"GROUND TRACK 

360 LOCATE 15,12 

370 PRINT"GEOSYNC 

poo LOCATE 17,12 

590 LOCATE 20,10 

400 PRINT"TO RUN THE PROGRAM 

PEO EOGALTE 21,10 

420 PRINT"TYPE FIRST TWO LETTERS 
eo 

caO! 

450 ‘accept choice from keyboard 
460 ' 


4] 


470 
480 
490 
500 
> 10 
520 
Do) 
540 
550 
560 
aeLO 
580 
D0 
600 
610 
620 
630 
640 
650 
660 
670 
680 
690 
700 
710 
720 
730 
740 
750 
760 
140 
780 
790 
800 
810 
820 
830 
840 
850 
860 
870 
880 
890 
900 
910 
920 
930 
940 
950 
960 
970) 
980 


X$=INPUT$(2) 
t 
‘load the chosen program 


IF X$="ic" THEN LOAD"INCL",R 

TF X$e ec THENILOnD 7ee an 

IF X$="1le" THEN LOAD"LEM",R 

IF X$=s"in" THEN LOAD ORE or 

IF X$="or" THEN SEOADU2BDORDI ak 

IF X$="gr" THEN LOAD"ELLIPORB",R 

IF X$a"ge" THEN LOAD"GEOSYNC",R ELSE 470 


END 

t 
fintroduct 1 on 
t 


EOGA Ei 23 ,0E 

PRD Nile These orbital programs demonstrate 
PRINT"several principals 6f orbitals tiligige 
PRINT 

PRINT"INITIAL CONDITIONS is a graphical 
PRINT"representation of the two-body 
PRINT"problenm. 

PRINS 

PRINT"ORBIT accepts the classical orbital 
PRINT"elements and draws the orbit around 
PRINT"a sphere representing the earth. 
PRINT 

PRINT"GROUND TRACK charts the ground trace 
PRINT"of the orbit over a map of the earth 
PRINT 

PRINT"GEOSYNC traces the ground track of a 
PRINT"geosynchronous orbit over a map of the 
PRINT"earth. 

LOCATE 24,1 

INPUT"press enter to continue....",KAL 

CLS 

RETURN 

SLOP 

t 


'title page 
t 


UNE 5 45 = oe Oa 

LOCATE 4,4:PRINT"SATELLITE ORBIT PROGRAMS eU LUZ nic 
LOCATE 6,6:PRINI"THE GRAPHICS CAPABILETES Cr Tie 
LOCATE 8,14:PRINT"MICROCOMPUTER 

LOCATE 11,19:PRINT"by 

LOCATE 14,10:PRINT"Kim Alldredge Langdorf 

LOCATE 16,15:PRINT"June 1986 

LOCATE 22,7 | 

INPUT"press enter to Containers ak A 

Gls 
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APPENDIX B 


INITIAL CONDITIONS PROGRAM LISTING 


10 ‘initial conditions 

20a: 

30°>"set up the sereemn 
q 


50 CLO:kK GY OFF: SGREE NM aye 
60 COLOR 9 
70 IF KAL<>0 THEN GOTO 110 
80 GOSUB 4810 

! 


100 ‘constants 
q 
120 PI=3.1415927# 


130 K=398603.2 
140 RE=6378.165 
LU 


0 

160 ‘counter to prevent redrawing of sphere 
170: 2 

ISO FIE A> 0 THEN | GOTOe 260 
190 ' 

200 'draw sphere 

210 F 

220 GOSUB 4450 

230 —" 

240 'set viewport on screen 
250 '! 


260. VY LEW C95, 30)-(310;5190) 72 
270 CLS:KAL=0 
! 


280 

290 ‘accept user input and echo on screen 

300 ' 

SVOMEOGATH 4, 1-PRINT™ ve RN a ":; LOCA 
7 ee een Ue eae elie ":LOCATE 10,1:PRae 
SP Rene See 18 12 PRN ";PR ONG 
SwbOGAGE. 21k PRIN To oe Ral Nels : 


32707 k0=03 VO= 0: BO=0:RAP= 0 

330° HOCATE 1.1 

240 INPUD"initial altitude in km ee 

350 RO=ALT+RE 

360 LOCATE “1,5- PRIN bi 
370 LOCATE 4,1:PRINT"Ro (km) = 

380° LOCATE 5,1:PRINT USING Fetter 

390 LOCATE 1,1 


400 INPUT"initial flight path angle in deg BE ADEG 


AVOVEGGATE2 ,VsPRINT” Me 
LOCATE 1,1 
420 IF BETADEG>90 GOTO 4120 
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Ao LOCATE 1,1:PRINT” 
1 


mie LOCATE 7,1:PRINI"Bo (deg) =" 

meme LOCATE 8,1:PRINT USING" #######.##' 5; BETADEG 
460 KAL=KAL+1 

470 BETA=BETADEG*(PI/180) 

wom, LOCATE 1,3 

490 PRINT"would you like to input initial 
mou LOCATE 2,3 

510 PRINT"speed, Vo, or desired apogee ,Ap? 
520 W$=INPUT$(1) 

530 IF W$="v" THEN GOTO 550 

meee iF WS$="a" THEN GOTO 4320 ELSE 520 


me 0 LOCATE els POR ole NOT 
".DRINT" MABOCATE 1,1 
bow) INPUT"initial speed in km/s = ",VO 


570 IF VO=0 THEN BEEP:GOTO 550 

poo LOCATE 10,1:PRINT"Vo (km/s)=":LOCATE 11,1:PRINT 
USING" #######.##" 5 VO 

pao LOCATE 1,1 

600 PRINT" 
10 

620 ‘calculations & type of orbit determination 

630 

640 VCIR=SQR(K/RO) 

650 VESC=SQR((2*K)/RO) 

660 SAM=VO*RO*COS(BETA) 

670 ENERGY=(V0*2/2)-(K/RO) 

680 ALPHA=(RO*V0%2)/K 

690 E=SQR((ALPHA—1)*2*(COS(BETA))*2+(SIN(BETA))“*2) 
mee Le ARPFHA=2 THEN SSS=3:GOTO 3160 

710 A=ABS(RO/(ALPHA-—2)) 

720 PP=A*(1-E%*2) 

730 IF ALPHA > 2 THEN SSS=4:GOTO 3090 

740 IF ABS(E-0)<.0001 THEN SSS=2:A=RO:R$="circular orbit": 
GOTO 790 

750 RPE=A*(1-E) : 

760 IF BETADEG=90 THEN SSS=6:GOTO 3240 

70) SSS=1:RAP=A*(1+E) 

moOe RS="elliptical orbit" 

790 T=SQR((4*PI"2*A%3)/K) 


800 

810 ‘to graphics routine 
820 ' 

830 GOTO 3650 

B40 S1OP 

850 ' 

860 ‘label strings 

S70. ' 

880 CLS 


890 A$="circular orbit" 
900 B$="elliptical orbit" 


oO 
9:20 
0 
940 
550 
960 
970 
980 
990 
1000 
IL e. 
LO20 
1030 
1040 
1050 
1060 
1070 


C$a" parabolic orbit, 
D$="hyperbolic orbit" 
E$s"circular orbit" 
F$="ballistic flight" 
G$="vertical launch” 
H$=s"semi-major axis 
I$="semi-latus rectum 
J$="eccentricity 
K$="alt at apogee 
L$="alt at perigee 
M$="vel at apogee 
N$="vel at perigee | 
O$="period 
P$="sp mech energy 
Q$="sp ang momentum 
S$="escape velocity 
T$="circular velocity 


-—-_ >» js. (j= 
—_ —_ 


i a nn an 
| —_— ll 


—_ 
> 


1080 
1090 


Us$="max alt above sfc 
V$="range 


——_—  —-_ 
—_> => 


1100 
PG, 
Lio 
1130 
1140 
1150 
1160 
1170 
Lis 
1190 
1200 
1 230 
1220 
E230 
1240 
eZ 
260 
LZ 
1280 
1290 
1300 
PS 
13720 
1330 
1340 
L330 
1360 
1370 
1380 
ew) 
1400 
1410 
1420 


W$="flight time 
XX$="vertical max alt 
YY¥$="rise time (min) 
UU$="alt above sfc 
HH$="Semi-minor axis 
? 


—_ 
= 


"display data 
f 

GOTO 
GOTO 
GOTO 
GOTO 
GOTO 
GOTO 


er 

IF 

iets 

ne 

IF 

re 
' 


9S$S=1 
$$$=2 
95953 
SSS=4 
555=5 
5$S=6 


1250 
Pr 7aLO 
1960 
2180 
2440 
2640 


"display 
' 


B=SQR(A° 2*(1-E°2)) 

LOCATE 4,10:PRINT Ks 
LOCATE 8, le PRUNE sais 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 


Se 2c e Rk iene 
7 le Nea 


7327 RUN ome 
O52 PRUNT es 


9,325 R ING Kae 
LO; PR iis 


eee ORS 


11,32 7PR UND om: 


—_> —-> 
—_> > 


elliptical orbit data 


8,20:PRINT USING ######## ##" 5B 


7,20:PRINT USING" ######## HSA 


9,20:PRINT USING ########.##" 3 PP 


10,20:PRINT USING"#####.#4###";3E 


11,20:PRINT USING "########.##";RAP-RE 
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1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
f10 
iy 2 O 
ies 0 
1540 
i> 0 
1560 
10 
LSeh®) 
190 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
0 O 
7 10 


LOCATE 12,1:PRINT L$ 
LOCATE 12,20:PRINT USING" ########.##";RPE-RE 
LOCATE 12,32:PRINT"km" 

MOCATE 13,1:PRINT M$ 

ae = SAM/RAP 

MOCATE 13,20:;PRINT USING" #####. #eoHH" VAP 
BOGATE 13,32:PRINT"kn/s" 

VPE= SAM/RPE 

POCGCATE l45:PRINT N$ 

LOCATE 14,20:PRINT USING" #####.#####";VPE 
LOCATE 14,32:PRINT"km/s" 

POGALE 5,1. P RINE 205 

LOCATE 15,20:PRINT USING" #####.#####":T/60 
POGATE 15,32:PRINT'min" 

WOGATE 16,1:PRINT P$ 

LOCATE 16,20:PRINT USING"#####4.#####"; ENERGY 
LOCATE 16,32:PRINT"km"CHR$(253)"/s"CHR$(253) 
mOCALTE 17,12PRIN?T Q$ 

LOCATE 17,20:PRINT USING"########.##";:SAM 
MOGATE 17,32:PRINI"km"CHR$(253)"/3s" 

POCATE 18,1:PRENT SS 

LOCATE 18,20:PRINT USING "##### .#####";:VESC 
LOCATE 18,32:PRINT"km/s 

BPOGATE 19,1:PRINI T$ 

LOCATE 19,20:PRINT USING "#####.#####":VCIR 
LOCATE 19,32:PRINT"km/s 


iL 7 70 


30 
1740 
1750 
1760 
770 
1780 
1790 
1800 


GOTO 2740 

t 

‘display circulammorbit data 
LOCATE 11,10:PRINT A$ 


OGATS@?S. 1:PRINT J$ 

LOCATE 13,20:PRINT USING"####.###8#"3E 
LOCATE 14,1:PRINT UU$ 

LOCATE 14,20:PRINT USING"######.##";RO-RE 
LOCATE 14,32:PRINT"km 

LOCATE 15,1:PRINT O$ 

LOCATE 15,20:PRINT USING "#####4.##48 #8": 


SQRC(4¥*PI*2¥*A*%3)/K) /60 


0 
1820 
ite 10 
1840 
1850 
1860 
1870 
1880 
1890 
oO 
1910 
20 
Fon 


EQCATE lSps2sPRINT min" 

LOCATE 16,1:PRINT P$ 

LOCATE 16,20:PRINT USING" #####.#####";-K/(2*A) 
EOG@wme loos 2ePRINT km CHRS$(253)"/s"CHR$(253) 
BOGATE -17,1:FPRINT.AO$ 

LOCATE 17,20:PRINT USING" ### #####.##";RO*SQR(K/RO) 
EOGAMEmI7 ,o2:PRINT km'CHR$(253)"/s" 

LOGATE  48,1:PRINT S$ 

LOCATE 18,20:PRINT USING" #####.#####";SQR((2¥*K)/RO) 
LOCATE 18,32:PRINT"km/s 

LOCATE 19,1:PRINT T$ 

LOCATE 19,20:PRINT USING" #####.### ##";SQR(K/RO) 
LOCATE 19,32:PRINT"km/s 


1940 
1950 
1960 
1970 
1980 
1990 
2000 
2010 
ZZ 0 
2030 
2040 
ZOO 
2060 
2070 
2080 
2090 
2100 
2 EEO 
2120 
2130 
2140 
yi 5). 
2160 
Z ie 0 
2180 
2190 
2200 
Zao 
2220 
2230 
2240 
L256 
2260 
2210 
2280 
2290 
2300 
Za@ 
2320 
2350 
2340 
2350 
2360 
230 
2380 
Zoo) 
2400 
2410 
2420 
2430 
2440 
2450 


GOTO 2740 

t 

"display parabolic orbit data 
t 


LOCATE 10,10:PRINT C$ 

LOCATE 12,1:PRINT J$ 

LOCATE 12,20:PRINT USING"#####.#48 #4 #"SE 
LOCATE 13,1:PRINT L$ 

LOCATE 13,20:PRINT USING" ####4##4.4#"3ABS(RPE-RE) 
LOCATE 13,32:PRINT"km 

LOCATE 14,1:PRINT P$ 

LOCATE 14,20:PRINT USING" #####44#.44";ENERGY 
LOCATE 14,32:PRINT"km"CHR$(253)"/s"CHR$(253) 
LOCATE 16,1:PRINT S$ 

LOCATE 15,1:PRINT 

LOCATE 15,20:PRINT USING" ######4#.44":SAM 
LOCATE 15,32:PRINT"km"CHR$(253)"/s" 

LOCATE 16,20:PRINT USING" 4######4.##"3VESC 
LOCATE 16,32:PRINT"km 

LOCATE 17,1:PRINT I$ 

LOCATE 17,20:PRINT USING" ######4#.4#"5PP 
LOCATE 17,32:PRINT"km 

GOTO 2740 

? 


"display hyperbolic orbital data 
t 


LOCALE Oe O- eRe Tes 

LOCATE il 2, leeR iN Is 

LOCATE 12,20:PRINT USING"##### #####"°E 
LOCATE. 1] 3 aieee Roni a: S 

LOCATE 13,20:PRINT USING ########.##";:ABS(RPE-RE) 
LOCATE 13,3222 R02 km 

LOGATE 4,1: PRUNES 

LOCATE 14,20:PRINT USING" "########.##"; ENERGY 
LOCATE 14,32:PRINIT"km CHRS( 2530) > = (Chih se 
LOCATE 15,32:PRENE'km CHRS (2 5a007 cu 

LOCATE, 15231: PR ia 

LOCATE 15,20:PRINT USING" ########.##"3SAM 
LOGATE [6. i- PR iNT 3s 

LOCATE 16,20:PRINT USING" ########.##";VESC 
LOCATE 16,32:PRINT"km 

LOCATE 17,1:PRINT"hyperbolic exc sp = 
HES=SQR(VO0*2-VESC“%2) 

LOCATE 17,20:PRINT USING"########.##":HES 
LOCATE 175322. PR Unde am 

LOCATE 18,1:PRINT I$ 

LOCATE 18,20:PRINT USING" ########.##";:ABS(PP) 
LOCATE @Te532- PRL a tom 

GOTO 2740 

t 


'display suborbital data 
t 


Zeaao LOCATE 10,10:PRINT F$ 

Zoe) Ss LOCATE 12,1:PRINT U$ 

2480 LOCATE 12,20:PRINT USING"######4##.##";RAP-RE 
22908 LOCATE 12,32:PRINT"km 

2500 LOCATE 13,1:PRINT V$ 

Zero LOCATE 13,20:PRINT USING" ######F#. HH"; RANGE 
tee LOCATE 13,32: PRINT"kn 

Pore LOCATE 14,1:PRINT W$ - 

yee) LOCATE 14,20:PRINT USING" ##### #4444" 3 FLTTIME/60 
peo LOCATE 14,32:PRINT"min 

PeoOe LOCATE 15,1:PRINT“optimum Bo for Vo = 

2570 LOCATE 15,20:PRINT USING" #####.#####";BOPT*(180/PI) 
Zee) LOCATE 15,32:PRINT"deg 

e900 GOTO 2740 

ZowO STOP 

Zoro ' 

2620 ‘display vertical launch data 

Zoe ' 

Boe LOGATE 10,10:PRINT R$ 

MoesO LOCATE 12,1:PRINT XX$ 

2660 LOCATE 12,20:PRINT USING" ########.##";:RMAX-RE 
MeO) LOCATE 12,32:PRINT"km 

moe) LOCATE 13,1:PRINT YY$ 

2690 LOCATE 13,20:PRINT USING" ######## .##";TRISE 
meo LOCATE 13,32:PRINT"min 

ZC 

2720 ‘choice of another satellite 

70 

Zee LOCATE 1,12PRINT" 

SEOCATE 1.,1 

2750 PRINT"another satellite? 

2760 X$=INPUT$(1) 

eo ir X$="n' THEN LOAD"pick",R 

feo LF X$="y° THEN GOTO" 50 ELSE 2740 

2790 STOP 

2800 ' 

2810 ‘subroutine to calculate true anomaly 

2620) ' 

Zea X=(ROFVO>2) /K 

Zoq@e Y=(ALPHA*SIN( BETA )*COS( BETA) )/(X*COS(BETA)*2-1) 
2850 THETAO=ATNCY) 

2860 RETURN 

28m) | 

2880 ‘choice of viewing orbital data, another 
2890 ‘satellite, or return to menu 

2900 ! 

2910 LOCATE 1,1 

2920 PRINT" uf 
2930 HLOGATE 251 

29405 PRINT" u 
292 0N LOCATE 1, 1 

2960 PRINT"again from same altitude?" 


2970 


2-928 0 


ySalNPURSGy 


IF Ke Gers ya a an L ORCA ae 1 5. 1s: .PaReieN tae 


1GOL0 390 


Zoo 
3000 
3010 
3020 
3030 
3040 
3050 
3060 
S107 10 
3080 
3090 
3100 
Sao 
51720 
5 TSC 
3140 
3 5.0 
3160 
370 
3180 
SILC 
300 
3210 
3220 
Se 
3240 
5250 
3300 
S260 
3270 
5250 
3290 
5300 
ao 
33:20 
3.350 
3340 
S)0)5)8, 
3360 
33:70 
5380 
5370 
3400 
3410 
3420 
3430 
3440 
3450 
3460 


IF X$="n" THEN 3000 ELSE 2950 

LOCATE 1,1 | 
PRINT"would you like to see orbital data? 
Y$=INPUT$(1) 

IF Y$="y" THEN SCREEN 0,0:COLOR 9:GOTO 880 
IF Y¥$sa"n" THEN 27408 ELSE SsO0e 


SCE 

t 

‘calculations for hyperbolic orbit 
t 


R$="hyperbolic orbit" 
RPE=A*(E-1) 

GOTO 3650 

SLOE 

t 


‘calculations for parabolic orbit 
t 


R$="parabolic orbit" 
RPE=.5*RO*ALPHA¥(COS(BETA))°*2 
PP=RPE*(1+E) 

GOTO 3650 

> LOP 

t 


‘calculations for vertical launch 
t 


R$="vertical launch" 
IF ALPHA >= 2 THEN LOCATE 2,15:PRINT"vert escapes":GOTO 


RMAX=(2*RO)/(2-ALPHA) 

F=SQR(RO/RMAX) 

Q=ATN(F/SQR(-F*F+41)) 
TRISE=.5*SQR(RMAX"3/(2*K) )*(PI+SIN(2*Q)-2*Q) 


GOTO 3650 

SLOP 

t 

‘calculations for ballistic trajectory 
' 

R$="suborbital flight 

Soo=5 


RAP=A¥*(1+E) 
X=(1/E)*(ALPHA*COS(BETA) *2-1) 

IF X=l THEN X=.999999 

IF X=-1 THEN X=-.999999 
THETA1=-ATN(X/SQR(-X*X41)) 41.5708 
Y=(1/E)*((ALPHA*RO/RE*COS(BETA)*2)-1) 
IF Y=l THEN Y=.999999 

IF Y=-1 THEN Y=-.999999 
THETAE=-ATN(Y/SQR(-Y*Y41))+1.5708 
RANGE=(2 * PI - THETA1-THETAE) * RE 


50 


3470 
3480 
3490 
3500 
53) @) 
a 20 
PU 
3540 
a 20 
55.00 
Boy 0 
3 0 
5590 
3600 
3610 
3620 
3630 
3640 
3650 
3660 
3670 
3680 
3690 
3700 
S710 
oy 20 
6730 
3740 
57 50 
3760 
B77 0 
5730 
5790 
3800 
2310 
3820 
3830 
3840 
3850 
3860 
3870 
3880 
3890 
5200 
32306 
3920 
3930 
3940 
550 
3960 
aoe 
3980 


Z=(1/E)*(1-(RO/A)) 

IF Z=-1 THEN Z=-.999999 

IF Z=1 THEN Z=.999999 
EO=-ATN(Z/SQR(ABS(-Z*Z+1)))+1.5708 
Q=(1/E)*(1-(RE/A)) 

IF Q=-1 THEN Q=-.999999 

IF Q=1 THEN Q=.999999 
P=-ATN(Q/SQR(ABS(-Q¥*Q+1)))+1.5708 
EE=2*PI-P 
FLTTIME=SQR(A“*3/K)*(EE-E*SIN(EE)-EO+E*SIN(EO)) 
Pl=2/((1/(2*A-RE))+(1/(2¥*A-RO))) 
QQ=SQR(K*P1)/(RO*VO) 

IF QQ=-1 THEN QQ=-.999999 

IF QQ=1 THEN QQ=.999999 
BOPT=-ATN(QQ/SQR(ABS(-QQ*QQ+1)))+1.5708 
GOTO 2880 | 

STOP 


"graphics subroutine 
? 


PPP=2*PI+,.1 

hOGNEE 1,15 

PRINT R$ 

Mies, <> 1 THEN GOTO 3720 

GOSUB 4200 

GOSUB 2800 

SSS=6 THEN PPP=PI/2 

BETA>THETAO THEN THETAO=THETAO+PI 
BETADEG=90 THEN O000=PI/2:PPP=PI*1.5:GOTO 3800 
BETADEG<>0 AND ALT<>0 GOTO 3800 
BETADEG<>0 THEN THETAO=THETAO+PI 
RPE<RE-.1 THEN THETAO=THETAO+PI 

000=0 

FOR THETA=000 TO PPP STEP PI/30 

rane 

IF SSS=1 THEN R=(A*(1-E*2))/(1+E*COS(THETA)) 
iecoo=? THEN R=A 

IF S$SS=3 THEN R=(2*RPE)/(1+COS(THETA)) 

IF SSS=4 THEN R=(A*(E*2-1))/(1+E*COS(THETA)) 
IF SSS=6 THEN R=RMAX/RE:Y=0:GOTO 3910 

IF SSS=5 THEN R=A*(1-E*2)/(1+E*COS(THETA)) 
R=R/RE 

Y=R*SIN(THETA-THETAO) 

Y=Y*(215/160)*(5/6) 

X=R*COS(THETA-THETAO) 

IF SSS=3 AND Y>RD THEN GOTO 4090 

IF SSS=4 AND X<-RD THEN GOTO 4090 

IF SSS=3 AND X<-RD THEN GOTO 4090 

IF SSS=4 AND Y>RD THEN GOTO 4090 

'IF POINT(X,Y)=l1 THEN PRESET(X,Y),1:GOTO 10070 
PRESET(X,Y),1 

LINE-(X,Y),KAL+1,,&HAAAA 


39600 IF SSS=4250R  SS5=3 60 LOna0cd 

4000 PSET(X,Y),KAL+1 

4010 VEL=SQR(K*(2/(R*RE)-1/4)) 

4020 LOCATE 18,1:PRINT"Vo (km/s)= 

4030 LOCATE 19,1:PRINT USING" ######.###";3 VEL 
4040 LOGATE 21,1:PRIND aie (im 

4050 LOGATE 22, 1:PRINT USING" ###### FFF"; (R#*RE)- RE 
4060 IF KL=l1 GOTO 4080 

4070 IF R<799 THEN °SS5=52ECCA i a >: : PRINT" suborbital 
4080 NEXT THETA 

40907 1 EP Vp oo oeGe Loma so 

4100 GOSUB 2880 

Ase: 

4120 'limit Bo from 0 to 90 degrees 

4130. * 

4140 BEEP 

415) OC Aer 2 5 

4160 PRINT"choose an angle between O and 90 
4170 GOTO 390 

4 SOs LOOP 

4190 ' 

4200 'draw circle representing earth 

70 

4220 X0=0:YO=0:RD=(A/RE)*2 

4230 IF SSS=3 THEN RD=(PP/RE)*2 

4240 VIEW(95,30)-(310,190),,3 

4250 WINDOW(XO-RD, YO-RD)-(XO+RD,YO+RD) 

4260) CIRCLECO. GO) 4-os a1 

4270 PAINT(0,0),1 ~ 

4280 CIRGLE(RO/ RE, O)RRDy 4s kale 

4290 PAINT(RO/RE,O),KAL+1 

4300 Re URN 

453 80> > (OF 

4320 

4330 'calculate Vo from apogee height 

4340 ' 

43550 LOCATE (ike ener. 

PP Ree de To LOGA LE meet 
4360 INPUT “apogee altitude (km) = ",RAP 
4370 RAP=RAP + RE 

4380 X=RO/RAP 

4390 IF RO = RAP GOTO 2550 


4400 ALPHA=((2-X-X*2*COS(BETA)*2)-SQR((2-X-X*2 
*COS(BETA)*2)*2-4*(1-X)*(1-X*2*COS( BETA) “2))) 
/(1-X*2*COS(BETA)*2) 


4410 
4420 
4430 


4440 


4450 


4460 


4470 


VO=SQR((ALPHA*K)/RO) 
GOTO 580 
SLO: 


'draw demonstration picture 
q 


CLS 


oy 


4480 
4490 
4500 
4510 
4520 
4530 
4540 
4550 
4560 
4570 
4580 
4590 
4600 
4610 
4620 
4630 
4640 
4650 
4660 
4670 
4680 
4690 
4700 
4710 
4720 
4730 
4740 
4750 
4760 
4770 
4780 
4790 
4800 
4810 
4820 
4830 
4840 
4850 
4860 
4870 
4880 
4890 
4900 
4910 
4920 
4930 
4940 
4950 
4960 
4970 
4980 
4990 


GUREEEC1 60, 10) se 5 5...) 
BAGNT( 1 6OpsO>) 1 

Pot (loo, LOa).3 

DRAW"r80;u40;3 

PSET(240,105),2 

DRAW"e403153r5;d5; 

BOCATE 7,28:PRINT"lLocal 

LOGATE 8,27:PRINT"horizon 

FOGALE 15,25 

PRINT"Ro 

LOCATE 8,36:PRINT"Vo 

fmcLe( 240,105),2,2 

PAUNTC241,104),2 

PATINECZ39,106),2 
Geez 4O, 105),20,3,.27#PI,.5*PI 
PSET(248,93) 

DRAW"r5;u5 

MOGATER 11,32 

PRINT"Bo 

BOCATE 1,1 

PRINT"A1t = height above surface of earth 
PRINT"Re radius of earth = 6378.165 km 
PRINT"Ro Alt + Re = initial radius 
bOGA hs 13,28:PRINT"AlLt 

BOGE 13,24:PRINT'Re 


EOCGAGE 4,1 
Beet Boe= initial flicht path angle 
PRINT"Vo = initial speed 


EOCATE 23,1 
INPUT"press enter 
CLS 

RETURN 

? 


POMmmcomtanuess.. os (1 


' introduction 


t 

LOCATE 2,10 

PRINT"INITIAL CONDITIONS 

EOGATE 4, 1 

PRINT"This program demonstrates the two-body 
PRINT"problem. Enter the initial values for 
PRINT"altitude, flight path angle, and either 
PRINT"speed or desired apogee height. 

ER Ne 

PRINT"The program then calculates various 
PRINT"characteristics of the resulting orbit 
PRINT"and plots the orbit around the earth. 
PRINT"You then have the option of viewing 
PRINT"the orbital data or plotting another 
PRINT"orbit on top of the first. 

PRINT 

PRINT"I£ the program isn't working properly, 
PRINT"press the Ctrl key and the Break key 
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5000 
5010 
5020 
50a0 
5040 
5050 
5060 
507.0 
5080 
3090 


PRINT"together and then press the F2 key 
PRINT"to restart. If you would like to make 
PRINT"a copy of the screen on the printer 
PRINT’press the shift key and the PrtSc key 
PRINT"at the same time. 

LOGALE, 24.7 
INPUT"press enter to continue...",LLL 

CLs 

RETURN 

SLO 
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10 
20 
50 
40 
50 
60 
70 
80 
90 
100 
ie) 0 
120 
B30 
140 
lb iO 
160 
170 
180 


200 
210 
220 
753)0 
240 
Za 
260 
270 
280 
Zo 0 
300 
310 
a2 0 


be > 6378.165 km":GOTO 260 


63 0 
340 
a0) 
360 
oO 
380 
S70 
400 
410 
420 
430 
440 
450 


Wor bit 
f 


APPENDIX C 


ORBIT PROGRAM LISTING 


"set up screen 
t 


Mey OPF:CLS:SCREEN 1,0 


COLOR 9 
t 


"introduction 
t 


GOSUB 2300 
t 


"set viewport 


VIEW(60,20)-(260,180),,3 
t 


‘constants 


PI=3.141593 
190 TWOPIT=2*PI 
TORAD=1.745329E-02 
TODEG=57.2957795# 
K=7.436575E-02:DTHETA=.004369 


t 


"routine to space input 


t 


GOSUB 2100 
t 


"accept user input & echo on screen 


t 


INPUT"semi-major axis (in km)=",H 
ER=H/6378.165 
Rico LEN BEEPRCOCN TH 2,2:PRINI"semi-major axis must 


1 


"draw globe 
t 


GOSUB 1770 


LOCATE 


EOCATH 4,17 PRINTa 


PRINT" 
t 
LOCATE 
PRINT" 
LOCATE 
PRINT" 
LOCATE 


Zee RINT" 
km 


Or) > Eales 
koe NT 
lee eRe 
oP RENT! 
Pelee RINT” 


=" "PRINT USING #####.#"3H 


DO 


PRINT" '.PRINT" "! 
LOCATE 19 9) een DN es 4 


PRINT" “SP RUN y 

GOSUB 2100: INPUT eecentricme,— ae 

LOCATE 8,1:PRINT"e =":PRINT USING" .######"5E 

GOSUB 2100 
INPUT"inclination, (decj= 0 

LOCATE 11,1:PRINT 2 ="PRENT USING ti ee 

PRINT" deg 

GOSUB 2100 

t 

'if eccentricity = 0 then skip arg of per and ascwmede 


’ 


[fF B=0 THEN "GOTO e710 

! 

INPUT"argument of perigee (deg)=",W 

LOCATE 15,1:PRINT"w =":PRINT USING"###.###";3 W 

PRINT" deg 

GOSUB 2100 

INPUT"ascending node (deg)=",01 

LOCATE 19,1:PRINT CHR$(234)" =":PRINT USING"###4.###"501 
PRINT" deg 

t 


"turn the sphere for edge view 
t 


01=01+89 
GOSUB 2100 
! 


"change to radians 
! 


W=W*TORAD 
T=I*TORAD 
O=01*TORAD 

! 


'P Q R coordinate system 
! 


PX=COS(W)*COS(0)-SIN(W)*SIN(0)*COS(I) 
PY=COS(W)*SIN(O)+SIN(W)*COS(0)*COS(I) 
PZ=SIN(W)*SIN(I) 
QX=-SIN(W)*C0S(0)-COS(W)*SIN(0)*COS(T) 
QY=-SIN(W)*SIN(0)+COS(W)*COS(0)*COS(I) 
QZ=COS(W)*SIN(I) 


"mean motion 
t 


N=K/ER“1.5 
f 

‘period 

t 
P=TWOPI/N 


GOSUB 2100 
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980 ' 

990 'main program loop 

oOo ' 

ome, FOR TTT=0 TO (P+.1) STEP P/50 
no20 ' 

1030 ‘kepler's equation 

read ' 

1050 GOSUB 1640 

1060 ' | 

1070 'xw yw rw coordinates 

1080 ' 

1090 XW=ER*(COS(ENI1)-E) 

1100 YW=ER*SQR(1-E*2)*SINCEN1) 

1110 RW=ER*(1-E*COS(EN1)) 

mee) RX=RW*6378.165 

ines 

1140 ‘calculate velocity & N (rev/day) 
50 ' 

1160 VEL=SQR(398603.2*((2/RX)-(1/H))) 
meee NI=16.997*(1/ER)°1.5 

1180 ' 

1190 'print vel, N, period, rad on screen 
mZO00 ' 

1210 GOSUB 2140 

2 ed 

1230 ‘determine if orbit hits earth 
ezeO 

1250 IF RW <= 1 THEN LOCATE 3,1:PRINT" Satellite will hit 
Garth (suborbital)":LOCATE 2,1 

mZ60 ' 

Maw ‘x y z coordinates 


1290 X=XW*¥PX+YW*QX 

1300 Y=XW*PY+YW*QY 

1310 Z=XW*PZ+YW*QZ 

1320 R=SQR(X*2+Y*24+Z%2) 

imo) Z1=Z*( 20/19) 

ieeoe LOCATE 23,1:PRINT Y POINT(-X,Z1) 

legoo ' 

1360 ‘check to see if orbit is behind earth 
lee.) ' 

1390 IF Y<O AND POINT(-X,Z1)<>0 GOTO 1490 
ee @ * 

20) pieot, orbit 

lao ' 

IanO Ct RGLE(—X,Z1),RO/20,2 

Wa 50" PAINT(-X,Z1),2 

1460 ' 

1470 ‘end of main loop 

1480 
1490 NEXT TTT 
L500: 


ar, 


150 
1520 
la 0 
1540 


leo 
1560 
1570 
1580 
oo 
1600 


"pause and then give choices 
t 
FOR S=l TO LOOQO:NEXT S 
LOCATE lean 

? 


LOCATE 1,5:PRINT"another satellite at same altitude? 
X$=INPUT$(1) 


IF. X$="n" THEN LOAD oaek {R 
IF X$="y" THEN A=0:1=0:E=0:W=0:01=0:GOTO 400 ELSE 570 
END 


LG ee 


1620 
1630 
1640 
Lo 0 
1660 
1670 
1680 
1690 
1700 
LNG 
1720 
0 
1740 
LoD 


"Kepler's eqn subroutine 
' 


EEE=90 
EEERAD=EEE*(PI/180) 

M=N¥*TTT ey se 

FOR I=1 TO 100 

IF ABS(EN1L-EEERAD)<.000001 GOTO 1740 

EEERAD=EN1 ; —_ 
MN=EEERAD-(E*SIN(EEERAD)) ime a oe 
DM=1-(E*COS(EEERAD)) 
EN1=EEERAD+((M-MN)/DM) 
NEXT I 

RETURN 

STOP 


oO 


L770 
LEeo 
ioe 
1800 
Po 
LS 2.0 
1830 
1840 
lene 
1860 
1870 
1880 
1890 
1900 
P30 
120 
Se 
1940 
oe) 
1960 
1970 
1980 
1970 
2000 
2OAe 


"routine to draw globe 
t 


,0=0-YO=0-hp= ek 

R=] 

A=PI 

COSA=COS(A):SINA=SIN(A) 

CES 

WINDOWCXO-RD, YO-RD)—-(X0+RD,YO+RD) 
CIRGCLECOSO) ees Sl 

PAINT CO-O) I 

FOR S=0 TO PI STEP PI/6 

X=0 

Y=R*=COSA*COsSes) 

RO=Roat NCS) 
CIRCLE(X, Y), R03, P1, fWOPI ASsee 
NEXT S 

FOR I=0 TO 6 

FOR J=0 TO 12 

T=I*PI/6 

S=J*PI/12 

COST=COSC )s:SINT—SiNc.) 
COSS=COS(S):SINS=SIN(S) 
X=R*COST*¥SINS 
Y=(R*¥SINA*SINT*¥SINS+R*COSA*COSS) 
LINE -(X,Y),3 
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2012 0 
2030 
2040 
2050 
2060 
2070 
2080 
209 0 
2100 
ZO 
ZO 
Za) 
2140 
mt 50 
2160 
2170 
2130 
Zo} 
2200 
we O 
2220 
Zee 0 
2240 
2250 
Pei) 
tet O 
2280 
me 9 
2310 
51) 
Ze 2.0 
2530 
2340 
2350 
2360 
Ze 10 
2530 
2390 
2400 
2410 
2420 
2430 
2440 
2450 
2460 
2470 
2480 


Both, 1) 53 


NEXT J 
NEXT I 
RETURN 
SOR 

f 


Fiproutine tor print bad, vel, per, .& N on screen 
t 


LOCATE 
PRINT" 
LOCATE 

RETURN 
t 


LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 


SS 


fT 


Les, 


1,1:PRINT"Rad= 

1,5:PRINT USING" "###### ##"3RX 
1,15:PRINT"km 

1,21:PRINT"Vel= 

1,25:PRINT USING "###. #####";3 VEL 
1,35:PRINT "km/sec 
Zee REN per = 

2,5:PRINT USING" "###### ##"3P 
a> PRINT min 


EOCAIE = 2, 21°PRINT” N= 

POGAUe 2ezaocPRINT USING #7? #7 #2#"-N1 
LOCATE 2,34:PRINT"rev/day 

RETURN 

STOP 

t 

fintroduction 

t 


LOCATE 5,18 
Bien ORBIT 


LOCATE 
ERAN 
PRoeN T" 
PRINT" 
PRINT" 
Renn 
PRINT" 
PRINT" 
PRINT" 
LOCATE 


or 

This program allows you to enter 
five of the six classical orbital 
elements. The sixth, time, is 
represented by starting the orbit 
at the perigee point and ending 
one complete orbit later. You may 
then display another orbit of the 
same altitude over the first. 

Ie 


Pievieeo, > )—(315,150),3,B 
INPUT"press enter to continue...",KAL 


CLS 
RETURN 
= LOR 
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APPENDIX D 


GROUND TRACK PROGRAM LISTING 


10 ‘ground track 

26% 

30 'set up the screen 

40 ' 

SO KEY OFF 3CES= SCRE ENieieee 
60 COLOR 9 

70 ! 

SO intigoduwesLen 

90 '' 

100 GOSUB 2690 

LTO 

120 'draw the map of the earth 
ESOm 

140 GOSUB 2290 

150 '' 

160 ‘constants 

170 ' 


160 Pla 326 19S 
190 TWOPI=2*P I 
200 TORAD=1.745329E-02 
210 TODEG=57.2957795# 
220 K=7.436575E-02: DTHETA=.004369 
230 RE=6378.165 
t 


250 ‘change color of satellite tracks 
t 


270 KAL=KAL+1 

280 IF KAL=1 THEN 
290 IF KAL=2 THEN 
300 IF KAL=3 THEN 
310 IF KAL=4 THEN 
320 IF KAL=5 THEN 
330 IF KAL=6 THEN 
340 IF KAL=7 THEN 

t 


my CG) GG) Go) 
oouuw vw vd wv a 
NO Wr dD WF bo 


360 ‘input the orbital elements 
f 


JOO COC AE alec 

390 PRINT" ” 
400 LOCATE 1,1 

410 INPUT"semi-major axis (in km)=",A 

420 IF A<RE THEN BEEP:GOTO 380 

430 LOCATE 6,1:PRINT"a = 

440 PRINT USING"#####":A 


450 
460 
470 
480 
490 
500 
5) 1U6, 
yO 
530 
540 


560 
D0 
580 
590 
600 
610 


630 
640 
650 
660 
670 
680 
690 
700 
710 
720 
730 
740 
750 
760 
E70 
780 
0 
800 
810 
820 
830 
840 
850 
860 
870 
880 
890 
900 
oo 
920 
G50 
940 


PRINT" km 

A=A/RE 

MOCATE 1,1 
PRINT" a 
meGaAtE 1. 1 

nimi eccentricrty= , 6 

LOCATE 10,1:PRINT"e = 

PRINT USING" .##";:E 

BOCATE 1,1 

PRINT" ‘ 
LOCATE 1,1 

PVEWM’ inclination (deg)=",I 

EeCGwwee io, l-PRINT i = 

PRINT USING"###4";I 

PRINT"deg 

Beeg= 0 THEN O1=0:W=0:GOTO 750 

MOCATE 1,1 

PRINT" . 
IOCATE Isl 

INPUT"argument of perigee (deg)=",W 
LOCATE 17,1:PRINT"w = 

PRINT USING"###";W 

PRINT"deg 

FOGATE 1,1 

PRINT" 
EOCATE 1,1 

INPUT"ascending node (deg)=",01 

LOCATE 21,1:PRINT CHR$(234)" = 

PRINT USING"###";01 

PRINT"deg 

LOCATE 1,1 

PRINT" " 
LOCATE 1,1 

INPUT"number of orbits=",NR 

EOCATE SP, 1: PRINT” 
? 


‘change to radians 
’ 


W=W*TORAD 
L[=I*TORAD 

O=01*TORAD 
t 


"change to P Q R coordinate system 
t 


PX=COS(W) *COS(0)-SIN(W)*SIN(O)*COS(TI) 
PY=COS(W)*SIN(0)+SIN(W)*COS(0)*COS(I) 
PZ=SIN(W)*SIN(L) 

QX=-SIN(W) *COS(0)-COS(W)*SIN(O)*COS(TI) 
QY=-SIN(W)*SIN(0)+COS(W)*COS(0)*COS(L) 
QZ=COS(W)*SIN(T) 


oo 


960 ‘calculate mean motion and period 

970 

980 N=K/A*%1.5 

990 P=TWOPI/N 

1000 ' 

1010 'draw the specified number of ground tracks 
FOZO*? 

1630 FOR TTTs0 TO NR*(P4.1) STEEP P7560 

1040 ' 

1050 'goto Kepler's equation 

1060 ' 

1070 GOSUB 2160 

LOSor 7 

1090 ‘change to xw- yw rw coordinates 

1100 XW=A*(COS(EN1)-E) 

1110 YW=A*SQR(1-E*2)*SIN(CEN1) 

1120 RW=A*(1-E*COS(EN1)) 

Lie. = 

1140 'is the semi-major axis smaller than the earth's 
radius? 

LiSsor* 

1160 IF RW <= 1 THEN LOCATE 1,1:PRINT" Satellite wil sone 
earth (suborbital)":FOR CCC=1 TO 1000:NEKT CCC:GOTO 2040 
hiv. 

1180 ‘change to x y z coordinates 

Lio 

1200 K=XW*¥PX+YW*QX 

1210 Y=XW*¥PY+YW*QY 

1220 Z=XW*¥PZ+YW*QZ 

1230" 

1240 ‘calculate the radius of the orbit 

1250 ' 

1260. R=SORCM On Oa) 

L270) * 

1280 ‘calculate the longitude and screen x-coordinate 
E290 * 

1300 ALPHA=ATN(Y/X) 

1310 ALPHA=ALPHA-DTHETA*TTT+0 

1220" * 

1330 'keep the correct quadrant 

PO 

1350 IF X<O AND Y>O THEN ALPHA=ALPHA+PI 

1360 IF X<O AND Y<O THEN ALPHA=ALPHA+PI 

1370 IF X>O AND Y<O THEN ALPHA=ALPHA+TWOPI 

1380 IF X>O AND Y>O THEN ALPHA=ALPHA+TWOPI 

1390 IF X=0 AND Y<=O THEN ALPHA=PI/2 

1400 IF X=0 AND Y>O THEN ALPHA=PI¥*1.5 

Ae 

1420 'change to degrees 
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1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
ie 0 
2 0 
30) 
1540 
£550 
O00 
coy O 
1580 
1590 
1600 
1610 
1620 
iro3 0 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
ielO 
1720 
7 3 0 
1740 
50 
O10 
1770 
1780 
1790 
1800 
foro 
1820 
e310 
1840 
ie 0 
1860 
1870 
1880 
1890 
1900 
too 
1920 


’ 


LONG=ALPHA*TODEG 

t : 
"keep longitude between O & 360 degrees 
t 


IF LONG>360 THEN LONG=LONG-360 
IF LONG<O THEN LONG=LONG+360 
t 


Healewlate the latitude 
LT=Z/R 
t 


"keep the latitude between O & 90 degrees 
t 


Peet >l THEN LT=1 

IF LT<-1 THEN LT=-1 
LT=ATN(LT/SQR(-LT*LT+1) ) 
t 


change latitude to degrees 

LAT=LT*TODEG 

add the ascending node to the longitude 
LONG=LONG-O1 

size the latitude and longitude to the screen 


LNG=LONG*(30/32) 
LNG=LNG*(32/36) 

LNG=LNG+169 

IF LNG>319 THEN LNG=LNG-300 
IF LNG<19 THEN LNG=LNG+300 
t 


"plot the groundtrack on the map 
' 


GURELECENG -100-LAT) , 196 
PAINT(LNG,100-LAT) ,C 
f 


"display the time, latitude and longitude on the screen 
f 


LOCATE 1,1:PRINT"TIME= 
LOGATE 1,6 

PRINT USING" #####.##" ; TTT 
EOCATE wl, 15:PRENT' min 
LOCATE 1,19:PRINT"LAT= 
LOCATE 1,23 

PRINT USING"###.##";LAT 
LOCATE 1,30:PRINT"LNG= 
LOGNLEN 1. 35 

PRINT USING" ###.##";LONG 


Pay 


1380 


1 


1940 ‘end of loop 

1950 ' 

1960 NEXT TTT 

1970 ' 

1980 'delay the last display 

1990 '' 

2000 FOR CCC=1 TO 1000:NEXT CCC 

2000) 

2020 ‘give option to do another satellite 
2030) 

2040 LOCATE 1) 1. PikerN i 

Ye LOG AEs Haak 

2050 PRINT"another satellite? 

2060 FOR, CCG=)) 20 1st - Ni kiacee 

2070) LOCATE ls 7 Rl Nile 

ee BOC AGE ele 

2080 X$=INPUT$(1) 

2090 IF X$="y" THEN GOTO 2396 

2100 IF X$="n" THEN LOAD picks hE Lot ez Ome 
Zoe 

2120 END 

2 aio 7 

2140 'loop to solve Kepler's equation 
ZO 7 

2160 EEE=90 

2170 EEERAD=EEE*(PI/180) 

230" MEN ait 

2190 FOR I=1 TO 100 

2200 IF ABSCENI-EEERAD)<.000001 GOTO 2260 
ZO eA 

2220 MN=EEERAD-(E*SIN(EEERAD)) 

2230 DM=1-(E*COS(EEERAD) ) 

2240 ENI=EEERAD+( (M-MN)/DM) 

2200 ONE al 

2260 RETURN 

2270, ShGr 

2730. 

2290 ‘subroutine to load map of the earth 
2300. ' 

2310 DIM SCRN( 4048) ,ATOS(6):LA=0 

2320 DEF SEG:BLOAD"arEryscrns VARe LRee sO) 
2330 PUTSCRN=VARPTR(ATOS(O) ):LA=VARPTR(SCRN(O) ) 
2340 BLOAD"picture",LA 

2350 CALL PUTSCRN(SCRN(0O) ) 

2500 REE N 

237033 EOP 

2330) 

2390 erase old elements & redraw grid 


2400 


2410 
2420 
2430 
2440 
2450 
2460 
2470 
2480 
2490 
2500 
ZO 
2520 
2 3.0 
2540 
290 
2560 
Ze 70 
2580 
Zo 2 O 
2600 
2610 
2620 
2630 
2640 
Zo 50 
2660 
2670 
2680 
2690 
27900 
Z710 
Ze 2 0 
D0 
2740 
2750 
2760 
ya 70 
2780 
2790 
2800 
Ze) 
2820 
2830 


MOCATE oO, | 

PRINT" y 

PRINT" ui 

PRINT" s 

EXCALE 10), 1 

PRINT" _ 

PRINT" 

ERIN, - 

EGOCGATE, 13,1 

PRINT" : 

PRINT" iy 

PRINT" i 

BOCGATE 17,1 

Een 

PRINT" p 

PRINT” a 

MOCATE 21,1 

PRIND® 2 

PRINT" u 

PRINT” G 

mINECL9,10)=(19;,190),1 

ENE Cro, 40)=(25,40),1 
LINE(19,100)-(25,100),3 

EeENE C193930)-( 25, 130), 1 

PONE Ge, 160)-(25,160),1 

GOTO 270 

SOP 

t 

‘int rodiac t ion 

t 

BOCATE 5,13 

PRINT'GROUND TRACK 

E@CATE 8. 1 

PRINT" This program plots the ground track 
PRINT" of an orbit on a cartesian map of 
PRINT" the earth. ‘Each grid is 30 degrees 
PRINT" by 30 degrees. The equator and 
PRINT" prime meridan are outlined in 
PRINT" white. 

PRINT 

ite 3 5 5 )—C505,150) ,3,8B 

RETURN 

> )Or 


AS 


APPENDIX && 


GEOSYNC PROGRAM LISTING 


10 geosync 

30 'set up screen 

ae OES Bb OFE2SCGREEN Vike 
70 PI=3.1415927# 


90 ‘introduction and draw map 

160 * 

110 GOSUB 1120 

120 GOSUB 1020 

130°" 

140 ‘accept and echo to screen user input 

1 5:07 

160° LOGATE Aga 

170 INPUT"inclination angle in degrees = ",III 
180 IF ILI>90 OR IT1<0 THEN BEEP GORGs Eo 
1[9OmMLOGATE See RENT ra 
2008 LOGATE, 1a 

210 INPUT"longitude (0 to 360 degj= JEG 

220 IF LG>360 OR LG<O THEN BEEP:GOTO 200 

230 XXX=LG*(30/36) 

240" LOGATES 151 

250 PRINT"inclination="II["deg longitude="LG"deg 
260 III=III*PI/180 

210F | 

280 'calculations for rotation of the earth 
290) 

300 WE=2*PI/24 

310 WO=WE 

320". 

330 'main program loop 

Sen 

350) FOR TOTO: lenses 

360 Y=SINCIII)*SIN(WO*T) 

370°? 

380 ‘calculate latitude 

5910. 

400 IF Y=1 THEN Y=.999999 

410 IF Y=-1 THEN Y=-.999999 

420 LAT=ATNC(Y/SQR(-Y*Y+4+1)) 

1 One 

440 ‘calculate longitude 

LeOsr" 

460 LONG=ATN(COS(III)*TAN(WO*T))-WE*T 

470 IF LONG < -PI/2 THEN LONG = LONG + PI 
480 IF LONG < O AND LAT <= -III THEN GOTO 710 


490 
500 
p10 
Bye 
530 
540 
550 
560 
0 
580 
DNG 
600 
610 
620 
630 
640 
650 
660 
670 
680 
690 
700 
710 
iz 
ae 
740 
750 
760 
He? O 
780 
790 
800 
810 
820 
530 
840 
850 
860 
870 
880 
890 
900 
910 
y20 


' 


"change to degrees 
t 

LAT=LAT*180/PI 
LONG=LONG*¥180*PI 

f 


"size to and keep on screen 

? 

LONG=LONG*(1/12) 
LONG=LONG+1694+XXX 3 

IF LONG<19 THEN LONG=LONG+300 
IF LONG>319 THEN LONG=LONG-300 
f 


"plot ground track 
t 


CIRCLE(LONG, 100-LAT),1,2 
PAINT(CLONG, l100-LAT),2 
' 


"end of loop 

f 

NEXT T 

t 

"4th quadrant fix 
’ 


FOR T=18 TO 24 STEP .25 
M=SINCLID)IASINCWO*T) 

IF Y=1l THEN Y=.999999 

IF Y=-1 THEN Y=-.999999 
LAT=ATNCY/SQR(-Y*Y+1)) 
LONG=ATN(COS(III)*TANCWO*T))-WE*T 
LONG=LONG+(2*PTI) 
LAT=LAT*180/PI 
LONG=LONG¥180*PI 
LONG=LONG*(1/12) 
LONG=LONG+169+XXX 

IF LONG<19 THEN LONG=LONG+300 
IF LONG>319 THEN LONG=LONG-300 
CIRCLE(LONG, 100-LAT),1,2 
PAINT(LONG, 100-LAT),2 

NEXT i 


"choice of another satellite 
t 


930 LOCATE 1,1:PRINT" 
PLOCATE 1,1 


940 
950 


PRINT"another satellite? 
FOR CCC=1 TO 1500:NEXT CCC 


960 LOCATE 1,1:PRINT" 
BoLOGA LE, 14.1 


970 
980 
290 


X$=INPUT$(1) 
IF X$="y" THEN GOTO 140 
ieee aan SHEN LOAD pick’,R ELSE 940 
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1000 
1010 
10720 
1030 
1040 
1050 
1060 
One 
1080 
1090 
i) 
iO 
ATG, 
PaO 
1140 
Leo 
io: 
1170 
1180 
1190 
1200 
1210 
EZ 
A S0 
1240 
1250 
1260 
1270 


SCE 

t 

‘subroutine to load map of the earth 
t 


DIM SCRN(4048),ATOS(6):LA=0 | 

DEF SEG:BLOAD"arryscrn",VARPTR(ATOS(0O)) 
PUTSCRN=VARPTR(ATOS(0)):LA=VARPTR(SCRN(0)) 
BLOAD"picture",LA 

CALL PUTSCRN(SCRN(0)) 

RETURN 

5 FGP 


1 


"introduction 
t 


LOCATE sql? 
PRINT"GEOSYNC 


EOGATE Ss | 

PRINT" A geosynchronous orbit is a circular 
PRINT" @enmbit located abowtw5,0 car enenadeer 
PRINT" 19,300 nm, or 35,750 km above the 
PRINT’ earth. Its period 4s 23"hapesommam 
PRINT 

PRINT" If a geosynchronous orbit has an 
PRINT" inclination of O degrees, it is 


PRINT" considered to be geostationary. 
LINE (335 )=63:157. 15005348 

RETURN 

STOP 


APPENDIX F 
SUGGESTED INPUT EXAMPLES 
fomstautethe programs, insert the disk ingthe "A" or 
default drive and turn on the computer. The program INTRO 
fame «Start with the title page Raear one on the screen, 
Press the enter key to et the program descriptions and 
then press enter again to arrive at the menu. To run a 
program, type the first two letters of the program name on 
the keyboard and the program will begin. 
The following input is suggested to familiarize the 
user with the programs: 
1, INITIAL CONDITIONS 


When the program asks for the initial altitude, 
enter "0", 


When the program asks for the flight path angle, 
enter "O", 


When the program asks whether you want to enter 
the initial speed or apogee altitude, type "v" for 


speed and then enter "9", 


After the program plots the orbit and asks if you 

want another satellite, type "n" for no. The 

program will ask if you want to see orbital data. 
7 Off 


Type “y for yes. 


The program will then ask if you want to try 


another satellite. Type "n" for no and the 


program returns you to the menu. 
Zee Orne le 


When the program asks for the semi-major axis, 
enter "8000". 
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on 


4. 


The program then draws@tne Spherpesrepresentinc es tie 
earth... It will then ask) porneuitemec ce) tcnetee 
Enter  s00oG ier 


The inclination is then requested. Enter "45", 


Enter "30" when the argument of perigee is asked 
fe Glee 


The last input. is the ascending node. Enter "60". 


The program then will draw the orbit around the 
sphere. When finished, the option to draw another 
orbit over the first is given. Typeliten” fore 


and to return to the menu. 


GROUND TRACK 


This program has the same input as the previous 
program, Enter the same data to see what the 
ground track of the orbit looks like. 


Another interesting orbit can be plotted by 
entering the following data: 


semi-major axis —2Go 20 lcm 
eccentricity = 4/0 
inclination = 65 degrees 
argument of perigee = 2/0 degrees 
ascending node = 85 degrees 
number of orbits = 2 


GEOSYNC 


There are only two inputs for this program. 
Enter "45" for inclination and) OU) econ stone: t lace 
Another interesting orbit can be obtained by 


entering "90" for inclination and either "0" or 
"TSO" “for longa tices 
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